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Abstract 

We consider a dilute Bose gas confined by a harmonic potential. We de- 
fine an appropriate thermodynamic limit and analyze the properties of the 
phases and phase transition in this limit. Critical properties in the pres- 
ence of the potential are found to be different from, though simply related, to 
those in the usual translationally invariant case. We argue that the properties 
of magnetically trapped rubidium Q and sodium Q gases (in which Bose- 
Einstein condensation has been recently observed) are well approximated by 
our thermodynamic limit except in a narrow window of temperature around 
the critical temperature. We also consider the effect of the confining potential 
on the non-equilibrium dynamics following a rapid quench to the ordered side 
and give a scaling description of the late time universal dynamics. 
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Much excitement has been generated in the last few months by the observation of Bose- 
Einstein condensation (BEC) in magnetically trapped alkali atoms |1|]. Since the original 
observation of BEC in rubidium atoms, two other groups have reported evidence for the same 
in lithium and sodium [[J atoms. A key feature of these experiments is the presence of a 
confining harmonic potential that the atoms feel. In this paper, we study the effects of such 
a confining potential on the statistical mechanics of ideal and interacting Bose systems (we 
confine our attention to repulsive interactions - so our results are not directly applicable to 
the Lithium experiment). In particular, we study the effects of this confining potential on the 
equilibrium critical properties of the finite temperature phase transition at which the Bose 
condensate first appears. We show how to define a sensible thermodynamic limit in which 
the critical singularities are present; the experiments are, of course, in finite systems - this is 
accounted for by finite-size scaling crossover functions which smooth out the singularites, in 
a manner quite analogous to that for familiar phase transitions in translationally invariant 
systems placed in a finite box. We will find that the exponents of the critical singularities 
are related, but not identical, to those of translationally invariant systems. The finite-size 
scaling crossover functions are expected to be new and not simply related to known crossover 
functions. We also consider the non-equilibrium dynamics associated with the establishment 
of long-range order following a rapid quench to the superfluid phase. We argue that recent 



results |L2[ obtained for the case of a translationally invariant system placed in a finite box 
may be used to understand the effects of the confining potential on this dynamics. 

On general grounds, there is no true phase transition to a Bose condensed phase in 
the presence of a confining potential for a finite number of particles. Such a transition 
is expected to appear only in a suitable "thermodynamic" limit; we will argue here that 
the experimentally appropriate limit is one in which the frequency, u, of the confining 
potential goes to zero, and the number of particles, N, to infinity, while keeping Nu 3 , the 
interaction strength, and the temperature (T) fixed. We will show that the critical properties 
in this limit can be understood by an application of the "local density approximation" of 
Oliva @]. In this approximation, the system is viewed as a collection of homogenous semi- 



macroscopic blocks, each with it's own local chemical potential. Each of these blocks can 
be treated independently of the others, and have properties characteristic of large uniform 
systems. The properties of the system with the confining potential can then be related to 
the corresponding properties of the usual translationally invariant systems. Many aspects of 
our general discussion below can be checked explicitly in a Hartree-Fock calculation |J . We 
thus obtain a complete understanding of the properties of the system in the thermodynamic 
limit. (In earlier treatments j|,f|, this approach was used for an approximate treatment 
of the non-critical properties of the inhomogenous Bose gases, but the importance of the 
thermodynamic limit in justifying it was not noted). For the phase transition, we estimate 
by the usual Ginzburg criterion that the system crosses over to the critical non-ideal regime 
only at small deviations from criticality (|T — T c \/T c ~ 10 -3 for the existing experimental 
systems); some properties of the system in this regime are quite significantly modified by 
the presence of the potential. Further, we show that the experimental systems are quite 
accurately described by our thermodynamic limit except in a narrow window of temperatures 
near the critical temperature, where finite-size crossovers need to be considered. However, 
this window is somewhat bigger than the temperature range for non-ideal behaviour; thus 
the crossover to non-ideal critical behaviour will be complicated by the presence of finite-size 
effects. 

We will describe our results using the Hamiltonian (in second-quantized notation) 



H= d 3 x 



|^-|W| 2 + Qmc V + y 2 + A 2 , 2 ) - H 2 + ||# 



(1) 



where m is the mass of the bosons, the interaction strength u is related to the scattering 
length a by u = Airh 2 a/m, and A is the anisotropy in the harmonic potential. When uj is 
finite, the system is confined to a finite region and will not have a true phase transition to 
a phase with long-range order. We expect "thermodynamic" behaviour to emerge only in 
the limit in which uj is sent to zero and the number of particles to infinity. In the usual 
case of a system confined in a box, this is achieved by keeping the density constant while 
sending the box size to infinity. For the oscillator case, this suggests that we scale the 



number of particles N with the volume over which the system is confined. This can be 
estimated at high temperatures (where interactions and quantum effects may be neglected) 
to be ~ {ksT /mu) 2 ) 'i . Thus, we guess that the thermodynamic limit is defined by sending 
uj to zero while keeping Nuj 3 fixed. As shown below, this does turn out to be a perfectly 
sensible limit in which there is a true phase transition at a finite non-zero temperature and 
the free energy per particle is finite. 

First consider the system in its high temperature phase. As uj goes to zero, the potential 
varies over a macroscopic length scale ~ 1/uj. Divide the system up into blocks, such that the 
potential does not vary significantly across any block. Each of these blocks is macroscopic 
in size and uniform. The correlation length in any block is some microscopic number (much 
smaller than the size of the block) determined by the local chemical potential. The different 
blocks can therefore be treated independently of each other. Furthermore, the properties 
of any block are well approximated by the thermodynamic behaviour of the corresponding 
uniform system. Extensive quantities such as the free energy are then a sum over the free 
energies of the individual blocks. For the total free energy F we thus have 

F= F biock = J d 3 rF(r), (2) 

blocks 

where T{f) is the free energy density of a block of size d 3 r centered at the point labelled 
by r. Note that we have replaced the sum over blocks by an integral, as the free energy 
density varies slowly from one block to another. The only r and uj dependence of F{f) 
is through the local chemical potential = /i — \muj 2 {x 2 + y 2 + \ 2 z 2 ) and therefore 
J-{r) = f{ujx,ujy,uj\z). Thus 

where x = ujx, y = ujy and z = \uz. Note that the function f{r) has no uj dependence. We 
immediately see that the free energy per particle F/N is finite in the limit iV — > oo, uj — > 0, 
Nuj 3 fixed. Thus the properties of the high temperature phase are related trivially to the 
corresponding properties of the uniform system. 
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This approach will fail if the correlation length in some of the blocks becomes bigger 
than the block size. For a small but non-zero u, this would happen close enough to a 
critical point. The true thermodynamic behaviour is accessed by first sending uj to zero and 
then approaching the critical point. When the limits are taken in this order, the block sizes 
can always be taken larger than the correlation lengths and so the equations (Q) and ([|) 
remain valid all the way upto a critical point. It is now easy to see that there is indeed a 
phase transition at a finite non-zero value of the temperature. As we have already seen, at 
high enough temperature there is a well-defined disordered phase; now consider the system 
at low enough temperature - there will always be some blocks with density bigger than 
the threshold value required for condensation at that temperature. These blocks will then 
become superfluid. Thus below a certain non-zero finite critical temperature the disordered 
phase is unstable to superfluid ordering. 

The properties of the system in the critical regime can be obtained straightforwardly 
from equations like (0) and ([$]). The singular part of the free energy density at point r 
satisfies the hyperscaling relation J~(r) = C/(£(f)) 3 , where £(f) is the correlation length in 
the block centered at f and C is a constant. For /i(f) sufficiently close to fi c , the critical 
point of the uniform system, £(f) ~ (/i c — [i(f))~ u where v ~ 0.67 is the correlation length 
exponent of the three-dimensional XY universality class. The total free energy is given by 

The integral is upto a value of r which is of the order of, but smaller than /mu 2 , and 

the ellipses denote terms regular in // — \i c . Note the extra factor of | in the exponent. This 
is a consequence of the quadratic potential in the system and leads to a violation of naive 
hyperscaling. Now T — T c is analytic in \i c — \i and hence the specific heat exponent 
is 2 - 30 + ±) ~ -1.52. This very weak singularity is likely to be masked by analytic 
background terms even if one could access the critical region. 

The order parameter correlation function for two points separated by x, in a block 
labelled by r, is G(f;x) ~ exp(— x/l;(r))/x l+v for x ^> £(r). The singular part of the 
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order parameter susceptibility, Xi (which is of physical importance in magnetic systems) 
varies as (T — T c )i~ L, ^ 2 ~ ri \ This exponent is positive (approximately 0.18) implying that the 
susceptibility is finite at the transition (even though there is a divergent correlation length 
associated with the transition) in contrast to the uniform case. 

We now consider the ordered phase. This phase is characterised by a non-zero expec- 
tation value for the Bose field The magnitude of this order parameter will be spatially 
inhomogenous (due to the confining potential), while its phase will be constant. First, we 
argue that there is indeed a well defined ordered phase in the thermodynamic limit defined 
above. Again, we imagine that the entire system is divided into semi-macroscopic blocks as 
before. We expect that both the density and the magnitude of the order parameter have 
short range correlations and hence their values in each block are determined by the local 
chemical potential ||. This immediately implies that the total number of particles scales as 
1/u 3 . Similarly, we also expect the free energy density to be determined just by the local 
chemical potential, and hence the total free energy also scales as 1/uj 3 . Thus the free energy 
per particle is again finite in the thermodynamic limit. 

Since the local chemical potential is a maximum at the center of the trap, the magnitude 
of the order parameter will be the largest at the center and will decrease as one moves away 
from the center. It eventually becomes zero when the local chemical potential becomes 
smaller than fi c . The position of this edge, denoted by f c , is given by fi(f c ) = fi c . Near the 
edge, ~ (M^*) — A*c) ~ ( r c — r) 13 where (3 ~ 0.34 is the usual order parameter exponent 
for the uniform system. This is an important point as it may be possible to study critical 
properties of the system by studying it close to the edge of the condensate in the ordered 
phase. 

The fraction of particles in the k = mode, n(k = 0) (which is of direct experi- 
mental significance) is proportional to (u 3 J d 3 r\ (&(r)) |) 2 in the thermodynamic limit. As 
one approaches the critical point from the ordered phase, the spatial extent of the con- 
densed region shrinks to zero as (// — /i c )2 ~ (T c — T)*. The maximum value of the 
order parameter (i.e. at r = 0) vanishes as (fi — fi c )P ~ (T c — T) 13 . This implies that 
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n(k = 0) ~ (T c — T) 2 ^ + i^ ~ (T c — T) 3,68 . An explicit calculation of the order parameter 
profile at T = within a Hartree-Fock approximation, was performed in Ref. [[|. An ex- 
tension of such a calculation to the critical regime |§ provides an illustration of the general 
discussion above (see also Ref. fllCfl ). 

As before, the free energy is determined completely once the free energy density of the 
uniform system is known as a function of [l. Thus the low temperature specific heat is 
proportional to T 3 , like in the uniform case, though the prefactor will be different. 

The phase of the order parameter, of course, has long range correlations. Associated 
with slow variations of this phase, we have the usual sound wave. The character of this 
mode is not different from the uniform case as long as the wavelength is small compared to 
a typical block size. The sound wave velocity will be determined in the usual way by the 
local superfluid density and compressibility. This spatial variation will cause the wave to be 
"refracted" as it propogates. 

The critical behaviour described above will become observable only in a window around 
Tq which is given by the usual Ginzburg criterion. This predicts that the system crosses 
over to the critical regime when £ ~ Aj- 2 /a where A^ = (h 2 /2mkBT)^ is the thermal de 
Broglie wavelength and a is the scattering length. This corresponds to \(T — T c )/T c \ ~ 
(a/Ay) 2 ~ 10 -3 — 10~ 4 for the Rb |l[ and Na || experiments. Outside of this window in 
the high temperature phase, the behaviour of the system can be well approximated by the 
thermodynamic limit of the ideal Bose gas in a harmonic potential and the calculations of 
fTT| are expected to apply. 

A small but non-zero value of uj leads to deviations from the theory described above, 
which is valid in the thermodynamic limit. "Finite-size" corrections are most significant 
near the critical point where the non-zero uj rounds off any singular behaviour of physical 
quantities. A crude estimate of the width of this region where finite-size effects are important 
may be obtained as follows: Near the bottom of the well, the correlation length varies 



appreciably over a length r ~ J {ji c — jj) / {muj 2 ). Finite-size effects will be negligible so 
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long as this length is much larger than £(r = 0). In the high temperature phase, by 
using the ideal bose gas expression for £(r = 0), this can be converted to the condition 
(T — T c )/T c < Two I '(ksTc) ~ 10~ 2 . A similar estimate is expected to hold below T c as well. 
Note that this window is at least an order of magnitude bigger than the temperature range 
where the system crosses over to the critical regime. Thus "finite-size" effects will prohibit 
observation of the true critical behaviour in the current range of experimental parameters. 
Increasing the scattering length and/or decreasing the frequency of the trap will enhance 
the possibility of measuring critical properties. 

So far we have restricted ourselves to the equilibrium properties of the system. How- 
ever, an interesting question which may also be experimentally relevant [0 is: How does the 
condensate grow in time to its final equilibrium value after a rapid quench in temperature 
from above T c to below ? This non-equilibrium question has been addressed recently [12| 



for a translationally invariant dilute Bose gas. It was shown that the standard phenomenol- 
ogy of phase-ordering kinetics []nj (which predicts, at late times, the existence of a single 
time dependent length scale ~ £*) can be used to obtain a scaling form for the equal-time 
correlation function of the Boson field \1/ : 

G(r,r',t,L) = (V>*(r,^(r',t)> = |(^)| 2 F( ±) (5) 

where L is the linear size of the system, c is a scale factor that depends on the final temper- 
ature (or equivalently, the final chemical potential), (i/j) is the equilibrium order parameter 
at that temperature and |r — r'\ is much larger than all microscopic scales in the problem. 
The value of z was numerically estimated to be close to 1 and it was argued that z should 
be exactly equal to 1 JTJ 

We now consider the effect of the confining potential on this dynamics. Imagine, as 
before, splitting the system up into many blocks (of size I ~ l/k>) such that the potential 
does not vary significantly within any block. As long as t* < I, the equal time correlation 
function within each block scales as: 

G(r, r', t,u) = | (ip(Ruj)} \ 2 H (|r - r'\u, c'(Ruj)tuj z ) (6) 



where R is the coordinate of the centre of the block, c' is a scale factor that depends on the 
local chemical potential and the value of z is unchanged from the translationally invariant 
case. Motivated by this we can now make the following ansatz for the scaling form valid 
for arbitrarily large values of t and \r — r'\ (in particular, r and r' may belong to different 
blocks): 

G = P(ruj,r'uj,tuj z ) (7) 

where z is the same as before. Note that as t — > oo P must reduce to (i^(r))(tp*(r')) which 
is consistent as (ip(r)) depends on r only through /i(uir). 

The total number of particles in the k = mode is proportional to u 3 J d 3 r J d 3 r' P and 
thus satisfies the scaling form: 

n(k = 0,t) =u~ 3 f(tu} z ) (8) 

The fraction of particles in the k — mode is given as n(k = 0,t)/N = (I / (N uo 3 )) f (tuj z ) . 
The function / of course depends on the value of temperature we quench to and the value 
of the anisotropy A. We thus conclude that the time dependence of n(k = 0)/N for fixed T, 
iVu; 3 and A but different values of u should exhibit scaling collapse at late times. For small 
values of its argument f(x) ~ ii For small enough uj (i.e for a large enough system) it may 
be possible to observe a corresponding window in time in which the condensate fraction has 
power-law growth n(k = 0,t)/N ~ if. 

In summary, we have presented a rather complete understanding of the effects of the 
confining potential on the equilibrium properties of a dilute Bose gas. In particular, we have 
demonstrated that the critical properties of the system in a confining potential are different 
from but simply related to the uniform case. We have also provided a scaling description 
for the late time universal dynamics following a rapid quench to the ordered phase. 
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